Effective Mass of a Vortex in a Clean Superconductor 
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We calculate the effective mass of a single quantized vortex in the BCS superconductor at finite 
temperature. Using the effective action for a vortex, we arrive at the mass formula as the integral 
of the spectral function J{uj)/lo' a over frequency. The spectral function is given in terms of the 
transition elements of the gradient of the Hamiltonian between Bogoliubov-deGennes eigenstates. 
We show that core-core, and core-extended transitions yield the vortex mass, near T = 0, of order of 
electron mass displaced by the normal core. The extended-extended states contributions are linearly 
divergent with the low frequency cutoff ui c , in accordance with the Ohmic character of the spectral 
function at low energies. We argue that the mass and friction are closely related in this system and 
arise from the same mechanism - interaction with the surrounding fermionic degrees of freedom. 



The dynamics of an individual vortex or their array is 
an integral part of our understanding of the superfiuid 
matter. In superconductors there have been phenomeno- 
logical approaches to write down a classical equation of 
motion of a vortex Q). This "effective" equation will in 
general rely on the information of the mass, along with 
the longitudinal and transverse forces that act on the 
vortex. 

The vortex mass also enters in an important way in the 
tunneling process (||| ■ Estimates of the mass based on 
Ginzburg-Landau functional gives a mass of order mkf 
per unit vortex length while more recent work em- 
ploying microscopic BCS-type analysis predicts a much 
larger mass ~ mfcf(fc/£ ) 2 ^he authors of Refs. 

^,[1 focused on the core quasiparticle contribution to 
mass, and were limited to zero temperature. We use in 
this work the effective action approach ||,|7j to go beyond 
those limitations. 

A quantized vortex is a phase singularity accompa- 
nied by a change in the order parameter density near 
the center. A closer scrutiny using the Bogoliubov-de 
Gennes (BdG) equation reveals that there are discrete 
energy levels inside the vortex core ||, in addition to 
a continuous spectrum with energy \E\ larger than the 
gap whose wavefunction extends beyond the core region. 
Both states are perturbed once the vortex is set in mo- 
tion, by a perturbing Hamiltonian equal to the inner 
product of the fermion current with the vortex velocity. 
Virtual transitions are induced among core and extended 
states and lead to an energy shift proportional to the 
square of the vortex velocity. Dividing by the squared 
velocity gives a mass in agreement with those of Refs. 

Here we proceed via an effective action in the Matsub- 
ara representation. By integrating out the fermions from 
the BCS action and expanding the remaining degrees of 
freedom up to second order in the vortex coordinate one 
can obtain the effective action of a single vortex ||[7j . It 
is possible to derive the transverse as well as the longitu- 
dinal forces on a vortex from this action M . Likewise we 



derive below a formal expression of the vortex mass valid 
at finite temperature. It agrees with the Hamiltonian 
perturbation theory at T — 0. Since core and extended 
states have non-vanishing overlap in the core region one 
expects core-core, core-extended and extended-extended 
transitions to contribute to the mass unlike previous the- 
ories which only considered core-core processes. Using 
approximate forms for eigenstates we explicitly compute 
the transition elements and the mass. At low tempera- 
ture the mass is comparable to that of electrons displaced 
by the hollow core of the vortex. 

The effective action by itself is non-local in time, hence 
dissipative, and the concept of mass only arises as a re- 
sult of the local-time approximation. We have deduced 
both the mass and friction of the vortex from a single 
framework, with the relative significance of inertia and 
dissipation largely determined by the motion itself. In 
this work we will focus on the pure system with no dis- 
order. 

The effective action of a single vortex, characterized 
by the two-dimensional coordinate r v (r) 0| in the BCS 
superconductor, is given by 
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In the theory of dissipative quantum systems J(w) is 
known as the spectral function ||: 

J( W )=^*(H^-^l)l/(^)-/(^)l IHV„!Zo|/3}| 2 ■ 
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It is the gradient of the Hamiltonian operator that gov- 
erns the longitudinal dynamics of the vortex. The trans- 
verse motion of the vortex, which is of course crucial for 
the overall dynamics, is not explicitly considered in this 
work. The longitudinal action shown in Eq. (Q) is respon- 
sible for the mass and also dissipation of the vortex. The 
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states a, (3 are the eigenstates of the Nambu Hamiltonian 
Ho with the energy E a and Ep respectively. 

By assuming r v (r) — r v (r') ~ (t — T')r v (r), one can 
rewrite the action as 
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The vortex mass follows as 
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for a homogeneous medium, a result known in dissipative 
quantum mechanics ]To[ |. 

The eigenstates of the BdG equation is written in the 
general form || 



ip a (x,y,z) = ~jjf kz 



(6) 



in the presence of a gap function A(r) = A(r)e for a 
singly-quantized vortex centered at the origin. The radial 
functions are obtained by solving the coupled differential 
equation 

r 2 u" + ru'+[(kr+E)r 2 -(n-l/2f]u = r 2 A(r)v 

r 2 v" + rv' + [(k 2 -E)r 2 -(p + l/2) 2 ]v = -r 2 A(r)u, (7) 

with k 2 = k 2 — k 2 . We set 2m = 1 for convenience and 
restore it in the final expression of M v . 

In arriving at the mass formula we made a physically 
important assumption that the dynamics of a vortex is 
local in time. From Eqs. ([!]) and (|^) one finds that 
such assumption is valid if the vortex motion is slow, 
i.e. r v (T)~r v (T'), on the relevant time scale \r— r'\. For 
situations where this assumption is violated, one must 
work with the original, non-local action. 

The eigenstates of the BdG equations are either local- 
ized at the core (\E\ < Aqq), or extended (\E\ > Aqq), 
where Aoo is a gap at infinity. Accordingly the transi- 
tion elements can be classified as core- core, core- extended 
(extended- core) and extended- extended. We examine the 
three contributions individually. The core solutions were 
first obtained by Caroli et al. [p| and the extended states 
are discussed by Cleary and by Bardeen et al. Q 

Core- core: The core states are given by 



u(r) 
v(r) 
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The state is normalized over the radius £o- The coherence 
length £o is much larger than k~ x in extreme type-H su- 
perconductors considered here. Thus the derivatives of 



the radial functions are essentially those of the Bessel 
functions, allowing us to apply the usual Bessel identi- 
ties. We obtain 

- (a\V v H \(3) = ^(x T m^ p ±x, (9) 

where eo = Aoo/fcr^o is the core level spacing and 
E a = [i a t{). From this we obtain the mass: 
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The factor kfL counts the number of electrons along the 
length of the vortex L, given that the electron spacing is 
kj . Since Ef/e ~ (fc/Co) 2 , the mass we obtained above 
is comparable to the mass of the electrons displaced by 
a hollow cylinder of radius £o. The zero temperature 
limit of Eq. (10) has been obtained previously d||. The 
temperature dependence is through the difference of the 
Fermi distribution function at ±Aoo, which is nearly con- 
stant for temperatures small compared with T c . Close to 
T c , the gap itself is sensitive to temperature, and Eq. jl0| ) 
diverges as £ ~ {T c -T)- 1 / 2 . 

This somewhat unphysical divergence of mass hints to 
us that perhaps our formula should not be taken too se- 
riously near the critical temperature. One possible rea- 
son is the transition element, Eq. (^|), which is based 
on Bessel function eigenstates. For divering coherence 
length, it is doubtful that the expression (||) can be taken 
seriously all the way to r = £ ^k~ 1 . 

The spectral function for the core is a delta peak at 
frequency to — eo- The (unapproximated) action due to 
the core-core transitions is 



S eff = 
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When eo ~ {T c — T) 1 / 2 is small, the exponential factor 
e ~a \r~T | remams order unity for a wide range of |r— r'|. 
For a large time interval the vortex may have moved a 
significant distance and not necessarily in a rectilinear 
fashion. In this case we cannot make the approxima- 
tion leading to Eq. (|3|). It still remains true, however, 
that Eq. (^) governs the dynamics of a vortex. A mo- 
tion described by Eq. (^) does not in general conserve 
energy. Hence a fast motion whose characteristic time 
scale is much smaller than €q 1 causes core quasiparticle 
transitions which dissipate energy. 

Core-extended: We choose a for the extended states, 
and (3 for the core states. For the extended states, there 
are two independent solutions for a given energy, whose 
approximate forms are 1 1 
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u(r) \ 




x J i(fc+r), 
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where «(£), w(-B) ={(E± v /E^A 3 (r))/2£J} 1 / 2 and fc± = 
(fc 2 ± y /E 2 -A 2 (r)) 1 / 2 . The normalization is over a large 
cylinder of radius R. The negative energy solutions can 
be found from the substitution (— v*, u*). Employing Eq. 
(H) for the mass, 



K«|v t 
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together with the constraint that the overlap integral is 
significant only if ^E^-A 2 ^ < A^ Q. The effective 
mass becomes (x = /3Aoo/2) 



A/,'; 
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sinh[x(l — /it)] 



„„M) (a) ^_ 

D7T / \ e o/J-i 1 — A* cosh [xj cosh [x^ij 



(14) 



with ft « 2.4. 

Replacing the integral with a characteristic value at 
fi = immediately shows that the core-extended vor- 
tex mass has the temperature dependence ~ /(— Aoo) — 
/(Aqo), similar to the core-core mass. We conclude that 
is of the same order as tf™, with the similar tem- 
perature dependence. The extended-core transitions will 
give the same mass as M£ e . 

Extended- extended: The matrix element appearing in 
Eq. (Q) consists of two parts. One involves the effect of 
"bending" of the order parameter dA(r)/dr and is only 
significant in the core region, < r < £o- The other transi- 
tion comes from the phase gradient Ve which falls off 
as 1/r, a long-ranged effect. One should also distinguish 
the transitions between states which lie on the same side 
of the Fermi surface, from those where a transition across 
the energy gap is involved. Their contributions to mass 
turn out to be very different. 

We consider the E a Ep > case first. Our approxi- 
mate wavefunctions Eq. (|l^) and their negative energy 
counterparts are supposed to be valid for r not much 
bigger than £o ■ As r gets well outside the core the wave- 
function picks up an additional phase shift not shown 
in Eq. (12) ||ll]| . Directly substituting our wavefunc- 
tions in the phase-gradient part of the integral leads to 
very small values for E a rj Ep, where the wavefunc- 
tions are likely to have the largest overlap. Hence a 
more accurate estimation of the wavefunction will not 
greatly enhance the size of the phase-gradient integral. 
We will thus assume that the integral can be restricted 
to the core region, where one can put A(r) = Aoo^/Co 
Q. With this substitution one has -(a\W v H\/3) w 
2ttA oc ^q 1 (x =p iy) j ° ru a (p)Vp( a )dr. We need to sepa- 
rately consider four cases where a and j3 are one of the 



extended eigenstates, and sum the squared magnitude of 
the transition amplitude for each process. We have 



|(a|V„ff„|/?>| 2 « ^(£ofc r -KI) 
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provided (E^-A^) 1 / 2 + (Bp-A^) 1 / 2 is on the order 
of Aqo or less. Both states lie close to the energy gap, 
reflecting the fact that states far from the Fermi surface 
do not influence the superconducting phenomenon. 

We can compute the spectral function J(u>) accord- 
ing to Eq. (||). The sum over states is with respect to 
k z , [i-a = ^±1, and the radial wavenumbers for both 
states. The integration measure is 
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When to = 0, the integral is logarithmically divergent be- 
cause of (E 2 —A 2 KJ )~ 1 term. Finite to effectively serves as 
the lower cutoff, giving 



J(to) 



Lm 2 C A 2 o (A oo +cj) 
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In a previous calculation by two of the authors (Ao 
& Zhu), it was shown that a vortex experiences friction 
at non-zero temperature The friction coefficient r\ is 
obtained as the derivative of J (to) at zero frequency [0: 



iA c 



k f L 



PA^/2 



4vr J cosh" 2 [(3 A x /2] 



ln(^] (IS) 

tO c 



with some lower frequency cutoff to c inserted by hand. A 
constant friction coefficient is a consequence of the Ohmic 
character of the spectral function, J(ui) ~ to (up to a log- 
arithmic factor). In turn one finds from Eq. ([|) that the 
mass in this case is linearly divergent! A straightforward 
evaluation of the mass gives 



Mf, e = m 



k f L\ A 



/3Aoo/2 



2tt 2 / uj c cosh 2 \J3 Aoo /2] 
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As anticipated, the mass diverges linearly for a fixed tem- 
perature as the low frequency cutoff tends to zero. For 
a fixed cutoff to c one gets exponentially vanishing mass 
in the zero temperature limit, because only the gapped 
transitions across the Fermi surface contribute to J (to). 

It was already mentioned that the very concept of mass 
depends on the motion the object executes. Here the 
mass is manifestly cutoff-dependent. To understand the 
physical origin of the cutoff to c let us go back to the orig- 
inal action, Eq. (Q). One can divide the integral over to 
into two parts; < to < io c and to c < to < oo, where a given 
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motion of a vortex is considered slow / fast in comparison 
with io~ r . For the fast spectral mode in w c <u< oo one 
can make the approximation r\>(r)— 7\,(t') « (j—T rs )f v (T) 
while for the more slowly responding environment one 
needs to keep the action as it is. If lu c is much smaller 
than Aqq, J{uj) is nearly Ohmic for lu<lo c , and the action 
becomes 
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This is the standard Caldeira-Leggett model of a Brow- 
nian particle ||], where the mass and friction are given 
by Eqs. (|ll]) and The friction coefheient is related 

to the mass by M£ e = 2r\j-KU) c . The classical equation of 
motion becomes 



M% e v + r]v = r\ \ v 

\1TUJc 



(21) 



Since the scale uj c is itself set by the motion, it is rea- 
sonable to assume that |v| w |w c v|, indicating that the 
effects of inertia and drag are always of the same order 
regardless of the characteristic time scale of the motion. 

Transitions across the Fermi level obviously do not con- 
tribute to friction, since the energy difference of the ini- 
tial and final states is greater than 2Aoo. We also expect 
that the mass contribution from this process is finite. In 
estimating the overlap integral we ignored the integra- 
tion in the far region, £o<r<R, because the phase shifts 
for states E a and Eg, E a + Eg w 0, should generically 
be different. Integrating over the core region gives 



(a|V„flo| 



A 2 
R 2 



A 2 

E 2 
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provided fc r £o > |/x Q | « \^g\. It is similar in magnitude 
to Eq. (H) for the E a E 8 >0 case. However the energy 
denominator is now of order (E a — Eg) 3 w A^, and free 
of divergence. The mass is given by m(kfL)[f(—A 00 ) — 
/(Aqq)] up to a numerical factor of order unity. While 
the temperature dependence is similar to M^ c and M£ e 
numerically it is much smaller than both. 

A vortex is not by itself assigned an intrinsic mass. 
Rather it is a consequence of interaction with the sur- 
rounding electrons. The motion of a vortex causes transi- 
tions between various BdG eigenstates and accrues mass 
as a consequence. 

We showed that both and M^ K ' are compara- 
ble to the mass of electrons displaced by the hollow core 
in the T = limit, and has a temperature dependence 
~ Co( T )[/(- A oo)-/(A 00 )]. The transition between ex- 
tended states with small energy differences leads to a 



frequency-dependent mass and a friction at finite tem- 
perature. Both quantities vanish at zero temperature 
exponentially. Considering extended states alone, one 
concludes that the vortex motion is dissipative at finite 
T even in the absence of impurities. In the low tempera- 
ture range where M^ c and M^ e are large and the friction 
is small, the vortex motion is Newtonian (mass times ac- 
celeration equals force), with a large mass. Experimental 
indications may have already existed [H . If the charac- 
teristic time scale of the motion is much smaller than the 
inverse core level spacing, however, the core level transi- 
tions can lead to dissipation also. 
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